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We review recent results concerning collective excitations in neutron-rich systems and reactions 
between charge asymmetric systems at Fermi energies. 

Solving numerically self-consistent transport equations for neutrons and protons with specific 
initial conditions, we explore the structure of the different dipole vibrations in the ^'^^ Sn system and 
investigate their dependence on the symmetry energy. We evidence the existence of a distinctive 
collective mode, that can be associated with the Pygmy Dipole Resonance, with an energy well 
^SJ \ below the standard Giant Dipole Resonance and isoscalar-like character, i.e. very weakly dependent 

on the isovector part of the nuclear effective interaction. At variance, the corresponding strength is 
' rather sensitive to the behavior of the symmetry energy below saturation, which rules the number 

CNj I of excess neutrons in the nuclear surface. 

<-«, , In reactions between charge asymmetric systems at Fermi energies, we investigate the interplay 

^ ^ between dissipation mechanisms and isospin effects. Observables sensitive to the isospin dependent 

' part of nuclear interaction are discussed, providing information on the symmetry energy density 

dependence below saturation. 



I. INTRODUCTION 



The Equation of State (EOS) of nuclear matter is of crucial importance for the understanding of a large variety of 
phenomena in nuclear physics and astrophysics. Transient states of nuclear matter far from normal conditions can be 
created in terrestrial laboratories and many experimental and theoretical efforts have been devoted to the study of 
nuclear reactions, from low to intermediate energies, as a possible tool to learn about the behavior of nuclear matter 
and its EOS. In particular, the availability of exotic beams has opened the way to explore, in laboratory conditions, 
new aspects of nuclear structure and dynamics up to extreme ratios of neutron (N) to proton (Z) numbers. Over the 

■ past years, measurements of isoscalar collective vibrations, collective flows and meson production have contributed 
to constrain the EOS for symmetric matter for densities up to five time the saturation value However, the EOS 
of asymmetric matter has comparatively few experimental constraints: The isovector part of the nuclear effective 

■ interaction (Asy-EOS) and the corresponding symmetry energy are largely unknown far from normal density. 

^vq ' This information is essential also in the astrophysical context, for the understanding of the properties of compact 
T-H ■ objects such as neutron stars, whose crust behaves as low-density asymmetric nuclear matter [2j] and whose core may 
^ I touch extreme values of density and asymmetry. Moreover, the low-density behavior of the symmetry energy also 
. ^ . affects the structure of exotic nuclei and the appearance of new features involving the neutron skin fsj . 

■ Over the past years, several observables which are sensitive to the Asy-EOS and testable experimentally, have been 



H . suggested 

- - ■ In this contribution we will focus on some effects related to the low-density sector of the nuclear EOS. The first 
part will be devoted to the study of collective excitations in neutron-rich systems. New exotic collective modes 
show up when one moves away from the valley of stability @. In particular, we will discuss results concerning the 
appearance of an interesting exotic mode, the Pygmy Dipole Resonance (PDR), which was observed as an unusually 
large concentration of the dipole response at energies clearly below the values associated with the standard Giant 
Dipole Resonance (GDR). Then, in the second part, we will discuss dissipation and fragmentation mechanisms at 
Fermi energies, concentrating on isospin effects and on the sensitivity of detectable observables to the Asy-EOS. 



II. TRANSPORT THEORIES AND SYMMETRY ENERGY 



Nuclear collective motion and nuclear reactions are modeled by solving transport equations based on mean field 
theories, with correlations included via hard nucleon-nucleon elastic collisions and via stochastic forces, selfconsistently 
evaluated from the mean phase-space trajectory, see |^. Stochasticity is essential in order to get distributions as 
well as to allow for the growth of dynamical instabilities. 
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In the beam energy range up to a few hundred MeV/u, the appropriate tool is the so-called Boltzmann-Langevin 
equation (BLE) 9]: 

| = ^ + {/,i7} = /.oH[/] + 5/[/], (1) 

where f{r,p,t) is the one-body distribution function, the semi-classical analog of the Wigner transform of the one- 
body density, H{r, p, t) the mean field Hamiltonian, IcoU the two-body collision term incorporating the Fermi statistics 
of the particles, and 6I[f] the fiuctuating part of the collision integral. Here we follow the approximate treatment to 
the BLE introduced in ,10], the so-called Stochastic Mean Field (SMF) model, where fluctuations are injected just in 
coordinate space by agitating the density profile. 

The symmetry energy, E^ym, appears in the energy density e{p, Pi) = e{p) + pEsym/A [pij pY -|- 0{pi/p)'^ + .., 
expressed in terms of total (p = pp + pn) and isospin (pi = Pp — Pn) densities. Egym gets a kinetic contribution 
directly from basic Pauli correlations and a potential part, C{p), from the highly controversial isospin dependence of 
the effective interactions: 

-± = + ^^^^'^ " T + 

The sensitivity of the simulation results can be tested against different choices of the density dependence of the 
coefficient C{p). 



III. COLLECTIVE EXCITATIONS IN NEUTRON-RICH SYSTEMS 



One of the important tasks in many-body physics is to understand the emergence of the collective features as well as 
their structure in terms of the individual motion of the constituents. The experimental characterization and theoretical 
description of new exotic collective excitations is a challenge for modern nuclear physics. Recent experiments provided 
several evidences about their existence but the available information is still incomplete and their nature is a matter 
of debate. In particular, many efforts have been devoted to the study of the PDR, identified as an unusually large 
concentration of the dipole response at energies below the values corresponding to the GDR. The latter is one of the 
most prominent and robust collective motions, present in all nuclei, whose centroid position varies, for medium-heavy 
nuclei, as SO^-^/^Mey. From a comparison of the available data for stable and unstable Sn isotopes a correlation 
between the fraction of pygmy strength and isospin asymmetry was noted JJJ. In general the exhausted sum-rule 
increases with the proton-to-neutron asymmetry. This behavior was related to the symmetry energy properties below 
saturation and therefore connected to the size of the neutron skin 0, [l^, [l^ . 

In spite of the theoretical progress in the interpretation of this mode and new experimental information (l3 - [T7| , a 
number of critical questions concerning the nature of the PDR still remains. Our goal is to address the important 
issue related to the collective nature of the PDR in connection with the role of the symmetry energy. 

An accurate picture of the GDR in nuclei corresponds to an admixture of Goldhaber- Teller (G-T) and Stenweidel- 
Jensen (S-J) vibrations. The latter, in symmetric nuclear matter, is a volume type oscillation of the isovector density 
Pi = pn — Pp keeping the total density p = Pn + Pp constant [l^. A microscopic, self-consistent study of the collective 
features and of the role of the nuclear effective interaction upon the PDR can be performed within the Landau 
theory of Fermi liquids. This is based on two coupled Landau- Vlasov kinetic equations (see Eq.(l), neglecting the 
stochastic term) for neutron and proton one-body distribution functions fq{f,p, i) with q — n,p, and was applied quite 
successfully in describing various features of the GDR, including pre-equilibrium dipole excitation in fusion reactions 
p^. However, it should be noticed that within such a semi-classical description shell effects are absent, certainly 
important in shaping the fine structure of the dipole response [2Q] . By solving numerically the Vlasov equation in 
the absence of Coulomb interaction, Urban [2l| evidenced from the study of the total dipole moment D a collective 
response around 8.6 MeV which was identified as a pygmy mode. It was pointed out, from the properties of transition 
densities and velocities, that the PDR can be related to one of the low-lying modes associated with isoscalar toroidal 
excitations, providing indications about its isoscalar character. Here, considering in the transport simulations also 
the Coulomb interaction, we shall investigate in a complementary way the collective nature of PDR by studying 
the dynamics of the pygmy degree of freedom, Dy, that is usually associated with the neutron excess in the nuclear 
surface [25|. Moreover, we shall explore the isoscalar character of the mode by a comparative analysis employing three 
different density parametrizations of the symmetry energy. 
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asy-EoS 


Esym I ^ 


L(MeV) 


i?n (fm) 


7?p(fm) 


ARnpi fm) 


asysoft 


29.9 


14.4 


4.90 


4.65 


0.25 


asystiff 


28.3 


72.6 


4.95 


4.65 


0.30 


asysupstiff 


28.3 


96.6 


4.96 


4.65 


0.31 



TABLE I; The symmetry energy at saturation (in MeV), the slope parameters, neutron rms radius, protons rms radius, neutron 
skin thickness for the three Asy-EOS. 



A. Ingredients of the simulations 



which we consider a Skyrme-like (SKM*) parametrization Uq = A— + B(—)"-^^ + C{p)— — —Tq + -- 



We neglect the two-body colUsion effects and hence the main ingredient of the dynamics is the nuclear mean- field, for 

Po Po ^ Po^ 2 dp po 

where Tq = -|-1(— 1) for q = n{p) and po denotes the saturation density. The saturation properties of symmetric nuclear 
matter are reproduced with the values of the coefficients A = — 356.8Mel^, B = 303.9Afey, a = 1/6, leading to a 
compressibility modulus K — 201MeV. For the isovector sector we employ three different parameterizations of C{p) 
with the density: the asysoft, the asystiff and asysuperstiff respectively, see Q for a detailed description. The value of 



the symmetry energy, Esym/ A, at saturation, as well as the slope parameter, L = 3po ^— — \p=po7 are reported in 

up 

Table U for each of these Asy-EOS. Just below the saturation density the asysoft mean field has a weak variation with 
density while the asysuperstiff shows a rapid decrease. Then, due to surface contributions to the collective oscillations, 
we expect to see some differences in the energy position of the dipole response of the system. 

The numerical procedure to integrate the transport equations is based on the test-particle (t.p.) method. For a 
good spanning of phase-space we work with 1200 t.p. per nucleon. We consider the neutron rich nucleus ^^^Sn and 
we determine its ground state configuration as the equilibrium (static) solution of Eq.(l). Then proton and neutron 

densities pq{r,t) — / - — 7-r^fq{^,p,t) can be evaluated. As an additional check of our initialization procedure, the 
J yZTTii) 

neutron and proton mean square radii (j'^) = J pq{r, t)(fir, as well as the skin thickness AR„p = V {r^ ) — ^J^J^^ 

were also calculated in the ground state and shown in Table HI The values obtained with our semi-classical approach are 
in a reasonable agreement with those reported by employing other models for similar interactions [23{ . The neutron 
skin thickness is increasing with the slope parameter, as expected from a faster reduction of the symmetry term on 
the surface [1, [l3| • This feature has been discussed in detail in Q . 

To inquire on the collective properties of the pygmy dipole we excite the nuclear system at the initial time t ~ — 
30 fm/c by boosting along the z direction all excess neutrons (Ne — 32) and in opposite direction all core nucleons, 
while keeping the 'cm.' of the nucleus at rest (Pygmy-like initial conditions). The excess neutrons were identified 
as the most distant = 32 neutrons from the nucleus CM. Then the system is left to evolve and the evolution of 
the collective coordinates Y, Xc and X, associated with the different isovector dipole modes (pygmy, core and total 
dipole) is followed for 600/to/c by solving numerically the equations (1). During the time evolution the number of t. 
p. escaping from the system corresponds, on average, to less than half a neutron, while the total energy conservation 
is satisfied within 1.5%. 



B. Results for dipole oscillations 

As shown in figure [1] apart from the quite undamped oscillations of the Y coordinate, we also remark that the 
core does not remain inert. We plot the time evolution of the dipole Dy, of the total dipole D and core dipole Dc 
moments, for two Asy-EoS. As observed, while Dy approaches its maximum value, an oscillatory motion of the dipole 
Dc initiates and this response is symmetry energy dependent: the larger is the slope parameter L, the more delayed is 
the isovector core reaction. This can be explained in terms of low-density (surface) contributions to the vibration and 
therefore of the density behavior of the symmetry energy below normal density: a larger L corresponds to a larger 
neutron presence in the surface and so to a smaller coupling to the core protons. We see that the total dipole D{t) 
is strongly affected by the presence of isovector core oscillations, mostly related to the isovector part of the effective 
interaction. Indeed, D{t) gets a higher oscillation frequency with respect to Dy, sensitive to the Asy-EOS. The fastest 
vibrations are observed in the asysoft case, which gives the largest value of the symmetry energy below saturation. 
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FIG. 1: (Color online) The time evolution of the total dipole D (a), of the dipole Dy (b) and of core dipole Dc (c), for asysoft 
(the blue (solid) lines) and asysuperstiff (the red (dashed) lines) EOS. Pygmy-like initial excitation. 



In correspondence the frequency of the pygmy mode seems to be not much affected by the trend of the symmetry 
energy below saturation, see also next figure [51 clearly showing the different nature, isoscalar-like, of this oscillation. 



For each case we calculate the power spectrum of Dy-. \Dy{Lo)Y 



Dy{t)e~"^*dt\'^ and similarly for D. The 



results are shown in figure [5J The position of the centroid corresponding to the GDR shifts toward larger values when 
we move from asysuperstiff (largest slope parameter L) to asysoft EOS. This evidences the importance of the volume, 
S-J component of the GDR in ^^^Sn. The energy centroid associated with the PDR is situated below the GDR peak, 
at around 8.5MeV, quite insensitive to the Asy-EOS, pointing to an isoscalar-like nature of this mode. Hence the 
structure of the dipole response can be explained in terms of the development of isoscalar-like (PDR) and isovector-like 
(GDR) modes, as observed in asymmetric systems Both modes are excited in the considered pygmy- like initial 
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FIG. 2: (Color online) The power spec- 
trum of total dipole (a) and of the dipole 
Dy (b) (in fm'^/c^), for asysoft (the blue 
(solid) lines), asystiff (the black (dot- 
dashed) lines) and asysuperstiff (the red 
(dashed) lines) EOS. Pygmy-like initial 
conditions. 



conditions. Looking at the total dipole mode direction, that is close to the isovector-like normal mode, one observes 
a quite large contribution in the GDR region. On the other hand, although the pygmy mode has a more complicated 
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FIG. 3: (Color online) The same as in fig- 
ure[2]but for a GDR-like initial excitation. 



structure 21], the Y direction appears closely related to it. Indeed a larger response amplitude is detected in the 
pygmy region, see figure 2 (bottom). 

To check the influence of the initial conditions on the dipole response, let us consider the case of a GDR-like 
excitation, corresponding to a boost of all neutrons against all protons, keeping the CM at rest. Now the initial 
excitation favors the isovector-like mode and even in the Y direction we observe a sizable contribution in the GDR 
region, see the Fourier spectrum of Dy in figure [H From this result it clearly emerges that a part of the Ne excess 
neutrons is involved in a GDR type motion and the relative weight depends on the symmetry energy: more neutrons 
are involved in the pygmy mode in the asysuperstiff EOS case, in connection to the larger neutron skin size. We have 
also checked that, if the coordinate Y is constructed by taking the Ny most distant neutrons (with Ny < Ne), the 
relative weight increases in the FDR region. In any case, since part of the excess nucleons contributes to the GDR 
mode, a low EWSR value is expected in the FDR region. Indeed, in the Fourier power spectrum of D in figure [31 a 
weak response is seen at the pygmy frequency. In the case of the GDR-like initial excitation, i.e. boosting all neutrons 
against all protons, we can relate the strength function to Im{D{uj)) and then the corresponding cross section 
can be calculated. Our estimate of the integrated cross section over the FDR region represents 2.7% for asysoft, 4.4% 
for asystiff and 4.5% for asysuperstiff, out of the total cross section. Hence the EWSR exhausted by the FDR is 
proportional to the skin thickness, in agreement with the results of [26| . 



IV. ISOSPIN EQUILIBRATION AND FRAGMENTATION MECHANISMS AT FERMI ENERGIES 

In this energy range, reactions between charge asymmetric systems are charecterized by a direct isospin transport 
in binary events. This process also involves the low density neck region and is sensitive to the low density behavior of 
Esym, see Refs.jl^, nil and references therein. Moreover, it is now quite well established that the largest part of the 
reaction cross section for dissipative collisions at Fermi energies goes through the Neck Fragmentation channel, with 
intermediate mass fragments (IMF) directly produced in the interacting zone in semiperipheral collisions on short time 
scales [sdj . It is possible to predict interesting isospin transport effects also for this fragmentation mechanism since 
clusters are formed still in a dilute asymmetric matter but always in contact with the regions of the projectile-like 
and target-like remnants almost at normal densities. 

Results on these mechanisms, obtained with the SMF model, are discussed below. 



A. The isospin transport ratio 



In peripheral and semi-peripheral reactions, it is interesting to look at the asymmetries of the various parts of 
the interacting system in the exit channel: emitted particles, projectile-like (FLF) and target-like fragments (TLF), 
and in the case of ternary (or higher multiplicity) events, IMF's. In particular, one can study the so-called isospin 
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transport ratio, which is defined as 



R 



P,T 



{x^ - x^) 



(3) 



with x'"^ = ^{x^ +x^). Here, x is an isospin sensitive quantity that has to be investigated with respect to equihbration. 
We consider primarily the asymmetry (3 = 1= {N — Z)/A, but also other quantities, such as isoscaling coefficients, 
ratios of production of light fragments, etc, can be of interest @. The indices H and L refer to the symmetric reaction 
between the heavy (n-rich) and the light (rt-poor) systems, while M refers to the mixed reaction. P, T denote the 
rapidity region, in which this quantity is measured, in particular the PLF and TLF rapidity regions. Clearly, this 
ratio is ±1 in the projectile and target regions, respectively, for complete transparency, and oppositely for complete 
rebound, while it is zero for complete equilibration. 




FIG. 4: Isospin transport ratios as a 
function of relative energy loss. Upper 
panel: separately for stiff (solid) and soft 
(dashed) Asy-EOS, and for two parame- 
terizations of the isoscalar part of the in- 
teraction: MD (circles and squares) and 
MI (diamonds and triangles), in the pro- 
jectile region (full symbols) and the tar- 
get region (open symbols). Lower panel: 
quadratic fit to all points for the stiff 
(solid), resp. soft (dashed) Asy-EOS. 



In a simple model one can show that the isospin transport ratio mainly depends on two quantities: the strength 
of the symmetry energy and the interaction time between the two reaction partners. Let us take, for instance, the 
asymmetry /3 of the PLF (or TLF) as the quantity x. At a first order approximation, in the mixed reaction this 
quantity relaxes towards its complete equilibration value, Peq = {Ph + /3l)/2, as 



PP.T 



+ (/3 



H,L 



(4) 



where t is the time elapsed while the reaction partners stay in contact (interaction time) and the damping r is 
mainly connected to the strength of the symmetry energy [28|. Inserting this expression into Eq.([3]), one obtains 
R^pj. = ±e-*/^ for the PLF and TLF regions, respectively. Hence the isospin transport ratio can be considered as a 
good observable to trace back the strength of the symmetry energy from the reaction dynamics provided a suitable 
selection of the interaction time is performed. The centrality dependence of the isospin ratio, for Sn -I- Sn collisions 
at 35 and 50 MeV/u, has been investigated in experiments as well as in theory |27| - |2 9l|. and information about the 
stiffness of the symmetry energy has been extracted from the analysis presented in [29| , based on the ImQMD model. 

Here we investigate more in detail the relation between charge equilibration, interaction times and thermal equi- 
librium. Longer interaction times should be correlated to a larger dissipation. It is then natural to look at the 
correlation between the isospin transport ratio and the total kinetic energy loss. In this way one can also better 
disentangle dynamical effects of the isoscalar and isovector part of the EOS, see 



28| 



It is seen in figure |4] (top) that the curves for the asysoft EOS (dashed) are generally lower in the projectile region 
(and oppositely for the target region), i.e. show more equilibration, than those for the asystiff EOS, due to the 
higher value of the symmetry energy at low density. To emphasize this trend, all the values for the stiff (circles) and 
the soft (squares) Asy-EOS, corresponding to different impact parameters, beam energies and also to two possible 
parametrizations of the isoscalar part of the nuclear interaction (with and without momentum dependence, MD and 
MI), are collected together in the bottom part of the figure. One can see that all the points essentially follow a given 
line, depending only on the symmetry energy parameterization adopted. It is seen, that there is a systematic effect 
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of the symmetry energy of the order of about 20 percent, which should be measurable. Moreover, we notice that the 
quantity i? is a rapidly decreasing function of the degree of dissipation, Eioss , reached in the collision. This can be 
explained in terms of dissipation mechanisms maily due to mean-field effects, as predicted by the SMF model. Indeed, 
according to a mean- field picture, a significant degree of thermal equilibrium (i.e. a considerable Eioss) would imply a 
rather long contact time between the two reaction partners, thus certainly leading to isospin equilibration, that needs 
a shorter time scale to be reached. The correlation suggested in figure |4] should represent a general feature of isospin 
diffusion, as expected on the basis of dominant mean-field mechanisms, and it would be of great interest to verify it 
experimentally. 

B. Isospin dynamics in neck fragmentation at Fermi energies 

In presence of density gradients, as the ones occurring when a low-density neck region is formed between the two 
reaction partners in semi-peripheral collisions, the isospin transport is mainly ruled by the density derivative of the 
symmetry energy and so we expect a larger neutron flow to the neck clusters for a stiffer symmetry energy around 
saturation [5]. This mechanism leads to the neutron enrichment of the neck region (isospin migration). This is shown 
in figure [5] (left), where the asymmetry of the neck and PLF-TLF regions, obtained in neutron-rich reactions at 50 
MeV/u, are plotted for two Asy-EOS choices. 

From the experimental point of view, a new analysis has been recently published on Sn+Ni data at 35 MeV/u by 
the Chimera CoUab.fs^. A strong correlation between neutron enrichement and fragment alignement (when the short 
emission time selection is enforced) is seen, that points to a stiff behavior of the symmetry energy (L sa 70 MeV), for 
which a large neutron enrichment of neck fragments is seen (left). 



0.25 I ^ 1 : 1 ^ i ^ i r 




Snll2 + Snll2 Pin Snl24 + Snl24 Pin 

FIG. 5: Left panel: asymmetry of IMF's (circles) and PLF-TLF (squares), as a function of the system initial asymmetry, for 
two Asy-EOS choices: asystiff (full lines) and asysoft (dashed lines). Right panel: Ratio between the neck IMF and the PLF 
asymmetries, as a function of the system initial asymmetry. The bands indicate the uncertainty in the calculations. 

In order to build observables less affected by secondary decay effects, in fig. 5 (right) we consider the ratio of the 
asymmetries of the IMF's to those of the residues (/3res) for stiff and soft Asy-EOS. This quantity can be estimated 
on the basis of simple energy balance considerations. By imposing to get a maximum (negative) variation of Esym 
when transfering the neutron richness from PLF and TLF towards the neck region, one obtains: 

PiMF _ EsymiPR) 
l^res -^sym(p/) 

From this simple argument the ratio between the IMF and residue asymmetries should depend only on symmetry 
energy properties and, in particular, on the difference of the symmetry energy corresponding to the residue and neck 
densities {pij and p/), as appropriate for isospin migration. It should also be larger than one, more so for the asystiff 
than for the asysoft EOS. It is seen indeed in figure 5 (right part), that this ratio is nicely dependent on the Asy-EOS 
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FIG. 6: Let panel: 
the reaction ^^*Sn 



Average total charge per event, associated with IMF's, as a function of the reduced rapidity, obtained in 



Sn at 50 MeV/u. Results are shown for ImQMD calculations at b = 6 fm (thick line) and b = 8 fm 
(thin line) and for SMF calculations at b = 6 fm (full circles) and b = 8 fm (open circles). A soft interaction is considered for 
the symmetry energy. Right panel: N/Z of IMF's as a function of the reduced rapidity. Lines and symbols are like in the left 
panel. Results corresponding to a stiff asy-EOS are also shown for ImQMD (dashed line) and SMF (crosses), for b=6 fm. 



only (being larger in the asystiff case) and not on the system considered. If final asymmetries were affected in the 
same way by secondary evaporation in the case of neck and PLF fragments, then one could directly compare the 
results of figure 5 (right) to data. However, due to the different size and temperature of the neck region with respect 
to PLF or TLF sources, de-excitation effects should be carefully checked with the help of suitable decay codes. 



C. Comparison with the predictions of different transport codes 



A detailed investigation of isospin equilibration has been recently undertaken within transport codes based on 
the molecular dynamics (QMD) approach 29J. In comparison to the transport models considered before, mainly 
describing one-body effects (such as the SMF model, see Eq.(l)), such approaches, where nucleons are represented as 
individual wave packets of fixed compact shape (usually taken as gaussians) , may lead to approximate descriptions of 
mean-field effects. On the other hand, fluctuations and correlations should be well reproduced, especially in the exit 
channel of multifragmentation events. As shown in Ref . [29j . where charge equilibration is investigated for Sn + Sn 
reactions at 35 and 50 MeV/u, the ImQMD code predicts a quite different behavior with respect to SMF: the isospin 
transport ratio exhibits a rather fiat behavior as a function of the impact parameter. This seems to indicate that, 
even in the case of central collisions, the contact time between the two reaction partners remains rather short, the 
dissipation mechanisms being mostly due to many-body correlations rather than to mean-field effects. Thus the more 
explosive dynamics would lead to the lower degree of isospin equilibration observed. 

To examine more in detail the origin of the observed discrepancies, results concerning IMF {Z > 2) properties, 
obtained with the SMF and ImQMD codes, are compared in figure 6. In the left panel, the average total charge per 
event, associated with IMF's, is plotted as a function of the reduced rapidity, for the reaction ^^'^Sn + ^^"^Sn at 50 
MeV/u and impact parameters b = 6 and 8 fm. From this comparison it is clear that in ImQMD a larger number 
of light IMF's, distributed over all rapidity range between PLF and TLF, are produced. On the other hand, mostly 
binary or ternary events are observed in SMF, with light IMF's located very close to mid-rapidity. Then the different 
reaction dynamics predicted by the two codes may explain the different results seen for isospin equilibration especially 
in semi-peripheral and central reactions (b ~ 4-6 fm). The fast ImQMD fragmentation dynamics inhibits nucleon 
exchange and charge equilibration. On the other hand, in SMF dissipation is dominated by mean-field mechanisms, 
acting over longer time intervals and leading to stronger equilibration effects. Results on the neutron content of the 
neck region are illustrated in the right panel of figure 6, that shows the global N/Z of IMF's as a function of the reduced 
rapidity. As discussed above, SMF calculations clearly predict a larger N/Z for IMF's produced at mid-rapidity, with 
respect to PLF and TLF regions (isospin migration effect). The effect is particularly pronounced in the case of the 
asystiff parametrization. On the contrary, ImQMD calculations predict a minimum of the N/Z ratio at mid-rapidity. 
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The reasons of the these differences need to be further investigated. 



V. CONCLUSIONS 



We have reviewed some aspects of the rich phenomenology associated with nuclear reactions, from which interesting 
hints are emerging to constrain the nuclear EOS and, in particular, the largely debated density behavior of the 
symmetry energy. Information on the low density region can be accessed in reactions from low to intermediate 
energies, where collective excitations and fragmentation mechanisms are dominant. 

We have shown, within a semi-classical Landau- Vlasov approach, the existence, in neutron rich nuclei, of a collective 
pygmy dipole mode determined by the oscillations of some excess neutrons against the nuclear core. From the transport 
simulations the PDR energy centroid for ^^"^Sn appears around 8.5 MeV ^ rather insensitive to the density dependence 
of the symmetry energy and well below the GDR peak. This supports the isoscalar-like character of this collective 
motion. A complex pattern, involving the coupling of the neutron skin with the core dipole mode, is noticed. This 
reduces considerably the EWSR acquired by the PDR, our numerical estimate providing values well below 10%, but 
proportional to the symmetry energy slope parameter L, that affects the number of excess neutrons on the nuclear 
surface. 

Concerning nuclear reactions at Fermi energies, we have discussed some results on isospin sensitive observables. 
In particular, we have concentrated our analysis on the charge equilibration mechanism (and its relation to energy 
dissipation) and on the neutron-enrichment of the neck region in semi-peripheral reactions. From the study of the 
latter mechanism, for which new experimental evidences have recently appeared [s?], hints are emerging towards a 
stiff behavior of the symmetry energy around normal density. This is compatible with recent results from structure 
data, see for instance the review article [ssj . 

Finally, we have also considered results of different transport codes (see section 4.3), undertaking a comparison of 
ImQMD and SMF predictions for selected observables. Due to the different role of mean-field effects, vs. fluctuations 
and many-body correlations, in the two codes, the description of the reaction dynamics is model dependent (see also 
[SJ]). This affects also isospin sensitive observables. Thus the whole reaction path needs to be carefully checked 
against experimental data in order to get decisive conclusions about the density dependence of the symmetry energy. 
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